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Abstract 

We study the factorization of the hypergeometric-type difference equation of 
Nikiforov and Uvarov on nonuniform lattices. An explicit form of the raising and 
lowering operators is derived and some relevant examples are given. 

1 Introduction 

In this paper we will deal with the so-called factorization method (FM) of the hyper- 
geometric-type difference equations on nonuniform lattices. The FM was already 
used by Darboux [14] and Schrodinger [27, 28] to obtain the solutions of differential 
equations, and also by Infeld and Hull [17] for finding analytical solutions of certain 
classes of second order differential equations. Later on, Miller extended it to difference 
equations [18] and ^-differences -in the Hahn sense- [19]. For more recent works see 
e.g. [5, 6, 11, 29, 30] and references therein. 

The classical FM was based on the existence of a so-called raising and lowering opera- 
tors for the corresponding equation that allows to find the explicit solutions in a very 
easy way. Going further, Atakishiyev and coauthors [5, 9, 6] have found the dynamical 
symmetry algebra related with the FM and the differential or difference equations. Of 
special interest was the paper by Smirnov [15] in which the equivalence of the FM and 
the Nikiforov et all theory [25] was shown, furthermore this paper pointed out that 
the aforementioned equivalence remains valid also for the nonuniform lattices that was 
shown later on in [20, 29]. In particular, in [29] a detailed study of the FM and its 
equivalence with the Nikiforov et al. approach to difference equations [25] have been 
established. Also, in [12], a special nonuniform lattice was considered. In fact, in 
[12] the author constructed the FM for the Askey-Wilson polynomials using basically 
the difference equation for the polynomials. In the present paper we will continue the 
research of the nonuniform lattice case. In fact, following the idea by Bangerezako 
[12] for the Askey-Wilson polynomials and Lorente [22] for the classical continuous 
and discrete cases, we will obtain the FM for the general polynomial solutions of the 
hypergeometric difference equation on the general quadratic nonuniform lattice x(s) = 
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c\q s + C2q~ s + C3. We will use, as it is already suggested in [9, 15], not the polynomial 
solutions but the corresponding normalized functions which is more natural and useful. 
In such a way, the method proposed here is the generalization of [12] and [22] to the 
aforementioned nonuniform lattice. 

The structure of the paper is as follows. In Section 2 we present some well-known 
results on orthogonal polynomials on nonuniform lattices [10, 25, 26], in section 3 we 
introduce the normalized functions and obtain some of their properties such as the 
lowering and raising operator that allow us, in Section 4, to obtain the factorization 
for the second order difference equation satisfied by such functions. Finally, in Section 
5, some relevant examples are worked out. 



2 Some basic properties of the ^-polynomials 

Here, we will summarize some of the properties of the g-polynomials useful for the rest 
of the work. For further information see e.g. [25]. 

We will deal here with the second order difference equation of the hypergeometric type 



o-(s) 



A 



Ax(s - \) lVx(s) 



(1) 



a(s) = a(x(s)) - %f(x(s))Ax(s - \), r(s) = f (x(s)), 



where V/(s) = f(s) — f(s — 1) and A/(s) = f(s + 1) — /(s) denote the backward and 
forward finite difference derivatives, respectively, a(x(s)) and f(x(s)) are polynomials 
in x(s) of degree at most 2 and 1, respectively, and A is a constant. In the following, 
we will use the following notation for the coefficients in the power expansions in x(s) 
of a(s) and f(s) 



a 



a(s) = a[x(s)] = yx 2 (s) + a'(0)x(s) + a(0), f(s) = f[x(s)] = f'x(s) + f (0). (2) 

An important property of the above equation is that the k-order difference derivative 
of a solution y(s) of (1), defined by 



Vk{s) q = 



A 



A 



Jq Ain(s) Ax fc _ 2 (s) "' Ax(s) 
also satisfies a difference equation of the hypergeometric type 



y(s) = A^y(s), 



a(s) 



Ax k {s-\) LVx fc (s) _ 



+ T k (s) 



Ay k (s) q 
Ax k {s) 



+ HkVk{s)q = 0, 



(3) 



where x k {s) = x{s + |) and [25, page 62, Eq. (3.1.29)] 



a(s + k) - a(s) + r(s + k)Ax(s + k-\] 



k-i 



Axfc_i(s) 



f'-k 



m=0 



Ar m (s) 
Ax rn (s) ' 



(4) 



It is important to notice that the above difference equations have polynomial solutions 
of the hypergeometric type iff x(s) is a function of the form [10, 26] 



x(s) = Cl (q)q s + c 2 (q)q- s + c 3 (q) = Cl (q)[q s + q' 3 ^} + C 3 (q), 



(5) 
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where c\, C2, C3 and q^ = ^ are constants which, in general, depend on q [25, 26]. For 
the above lattice, a straightforward calculation shows that T k (s) is a polynomial of first 
degree in x k (s) of the form (see e.g. [10]) 



a" 



Tk(s) = f' k x k {s) + f fe (0), f£ = [2k] q — + a q {2k)f, 



r fc (0) 



C 3 CJ' 



(6) 



(2[fc], - [2*;],) + ?'(0)[£;] ? + <*/(<*, (fc) - a,(2fc)) + f(0)a q (k), 



where the q-numbers [k] q and a q (k) are defined by 



_1 > 

q 2 — q 2 



a ? (/c) 



,2 + g 2 



(7) 



and [n] q \ are the (/-factorials [n] g ! = [1] 9 [2] ? • • • [n] q . 

Both difference equations (1) and (3) can be rewritten in the symmetric form 



and 



Ax(s-i) _ 
A 



a{s)p{s 



Ax k (s - 5) 



Vy(£) 
Vx(s)_ 



+ X n p(s)y(s) = 0, 



+ PkPk{s)yk{s) = 0, 



where p(s) and Pfc(s) are the weight functions satisfying the Pearson-type difference 
equations 



A 



Ax{s - \) 



— [a(s)p(s)\ = t(s)p(s), 



A 



Ax k (s - o) 



TT [o-(s)Pfc(«)] = T k (s)p k (s), (8) 



respectively. In [25] it is shown that the polynomial solutions of (3) (and so the poly- 
nomial solutions of (1)) are determined by the q- analogue of the Rodrigues formula on 
the nonuniform lattices 



A 



Ax fc _i(s) Ax(s) 



— P n (x(8)) q EE A(*>P B (X(5)) S " 



where 



vi n) /( s ) 



V 



Vx k+ i{s) Vx k+2 {s) Vxn(s) 
k-l 



m s ) 



/(*)■ 



(9) 



n 



^- J] \a q {n + m-l)r' + [n + m-l] q ^\ (10) 

m=0 ' 



Thus [25, page 66, Eq. (3.2.19)] 



B 



Pn(x(s)) q = -n-VWpnis), VW EE 



V 



n 

where p n (s) = p(s + n) a(s + k) and 



(11) 



fe=i 



A™ = -[n] q <J a 9 (n - l)r' + [n - — 



(12) 
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In this paper we will deal with orthogonal g-polynomials and functions. It can be 
proven [25], by using the difference equation of hypergeometric-type (1), that if the 
boundary condition 

a(s)p(s)x k (s - -) = 0, Vfc > 0, (13) 

2 s=a,b 

holds, then the polynomials P n (s) q are orthogonal, i.e., 

6-1 1 

^2P n (x(s)) q P m (x(s)) q p(s)Ax(s - -) = 5 nm dl, s = a,a + l,...,b-l, (14) 

s=a 

where p(s) is a solution of the Pearson- type equation (8). In the special case of the 
linear exponential lattice x(s) = q s the above relation can be written in terms of the 



Jackson ^-integral (see e.g. [16, 21]) f(t)d q t, defined by 

fZ2 rz 2 pz\ 

\ f(t)d q t = / f(t)d q t - / f(t)d q t, 

JZ! JO JO 



'Zl 

where 



rz 

/ f(t)d q t = z{l-q)J2f(zq k )q k , < q < 1, 
Jo 



k=0 

as follows: 



I" P n {t) q P m {t) q uj{t)d q t = 5 nm q l l 2 d 2 nl t = q s , u(t) = uiiq 1 ) = p(t). (15) 
J q a 

Notice that the above boundary condition (13) is valid for k = 0. Moreover, if we 
assume that a is finite, then (13) is fulfilled at s = a providing that a(a) = [25, §3.3, 
page 70]. In the following we will assume that this condition holds. The squared norm 
in (14) is given by [25, Chapter 3, Section 3.7.2, pag. 104] 



b-n-l 



d 2 n = {-l) n A n , n B 2 n Pn(s)Ax n (s - 

s=a 

There is also a so-called continuous orthogonality. In fact, if there exist a contour Y 
such that 

/ A[p(z)a(z)x k (z --)}dz = 0, V/fc > 0, (16) 
Jr 2 

then [25] 

J P n (x(z)) q P m (x(z)) q p(z)Ax(z -^)dz = 0, n^m. 

A simple consequence of the orthogonality is the following three term recurrence rela- 
tion: 

x(s)P n (x(s)) q = a n P n+1 (x(s)) q + j3 n P n (x(s)) q + i n P n -i(x(s)) q , (17) 

where a n , (3 n and 7„ are constants. If P n (s) q = a n x n {s) + b n x n ~ 1 {s) + • • • , then using 
(17) we find 

a n Q b n b n+ \ a n -i d\ 

an = , Pn = , In = ^ " V 18 ) 

On+l On a n+l a n « n _i 

To obtain the explicit values of a n , (3 n we will use the following lemma, -interesting in 
its own right- that can be proven by induction: 
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Lemma 2.1 

A^x n (s) = 



n 



x n - K (s) + c 3 [n 



[n ~ ljgj 
'[n-k-l] q \ 



(n — fc) 



n 



[n-fc],!"** 
In the case k = n — 1, it becomes 

A^-^x^s) = [n],!x n _i(s) + c 3 [n - 1],! (n - [n] q 
Now, using the Rodrigues formula (9) for k = n — 1, 



„n— fc— 1 



(«) + ■ 



(19) 



B n V 



Pn-l(s) 



Pn-l(s) Vx n (s 



-PAS) 



as well as the identities p n (s) = p n -\{ s + l) cr ( s + 1)) x n( s ) = ^n-i( s + 5) an d the 
Pearson equation (8) for p n -i(s), we find 



Thus 



A^P n (x(s)) q = An^BnTn-^s). 
■^n,n—\Bnj' n —\ 



n-1 , 

#n J| I a q (n + k- l)f +[n + k- l] q — \ , 

b n _ [n] g f n _i(0) 



n 



and 



So 



a. 



'n-1 



B r , 



+ c 3(Ng ~n). 

a q (n- l)f' + [n-l] 9 ^ B n \ n [2n\ q [2n + 1], 



and 



# n+ i (a 9 (2n-l)f / + [2n-l] 9 ^)(a g (2n)f / + [2n] g ^) £ n +i [n], A 2n A 2n +i 

a. = ^#=^) _ ["+ynw + C3(Wg + ! _ [n + iy . 

^n— 1 

Using the Rodrigues formula the following difference-recurrent relation follows [1, 25] 



a{s) 



Vx(s) 



An 



\n\„T, 



q'n 



22 

T n {s)Pn(x(s))q ~ P n +l{x(s)) q 
-t>n+l 



where r n (s) is given by (6), where the identity f' n = — ^ 2ra+1 has been used. 

[2n + l\ q 

Then, using the explicit expression for the coefficient a n , we find 

a{s) Y^k = ^TM Pn{x{s))q _ ^P n+1 (x(s)) q . 



Vx(s) 



[2n], 



(20) 



This equation defines a raising operator in terms of the backward difference in the 
sense that we can obtain the polynomial P n +i of degree n + 1 from the lower degree 
polynomial P n . 

From the above equation and using the identity V = A — VA, the second order differ- 
ence equation and the three terms recurrence relation we find [1] lowering-type operator: 



[a{s) + r(s)Ax(s - l)]^^k = l2^Lp n _ l{x{s))q+ 



Ax(s) 
\ A 2n 



[2n], 



[2n] q 
(x(s) - P n ) 



(21) 



Pn(x(s)) q . 
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The most general polynomial solution of the g-hypergeometric equation (1) corresponds 
to the case 

4 4 

a(s) = AY\[s - Si] q = Cq~ 2s JJ(<? S - q Sl ), A, C,not vanishing constants (22) 

i=l i=l 

and has the form [26] 

g Sl+S2+M ) ^l+ S3 +M ) gSl+S4+M > 9 > 9 ) > ( 2 3) 

where D„ is a normalizing constant and the basic hypergeometric series p 9 are defined 
by [21] 

i /oi,...,a r N (oi;g)fc---(or;g)fc ^ fc \ ( nfc # (fc -i)] p ' r+1 
6 l9 ...,6 p ^(6i;g) fc ---(6 P ;g) fc («;g) fc ^ M J 

and 

fc-i 

(a;9) fc = n( 1 -«9 m ), (24) 

m=0 

is the (/-analogue of the Pochhammer symbol. Instances of such polynomials are the 
Askey-Wilson polynomials, the g-Racah polynomials and big g-Jacobi polynomials 
among others [21, 26]. 

3 The orthonormal functions on nonuniform lattices 

In this section we will introduce a set of orthonormal functions which are orthogonal 
with respect to the unit weight [9, 15] 

Ms) = ^/MKPn(x{s)) q , (25) 

e.g. for the case of discrete orthogonality we have 



b-i 1 

^2 L Pn{ s i) ( fm{Si)^x{s i - -) = <W 



2' 

Si=a 

Next, we will establish several important properties of such functions which gener- 
alize, to the nonuniform lattices, the ones given in [22]. In the following we will use the 
notation 8(s) = a(s) + t(s)Ax(s — ^). 

First of all, inserting (25) into (1), (17), (20), (21) we obtain that they satisfy the 
following difference equation: 

Ve( s )<r( s + i)^L^ n ( s + 1) + ye( s -i))a( s )-l-^( s - 1)- 

(26) 

( @(s) a(s) \ . , , A . , , K ' 

the three term recurrence relation: 

an^p-(Pn+l(s) + Jn^f^Vn-l(s) + (/3 n - x(s))<p n (s) = 0, (27) 
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the raising-type formula: 



T + ( \ 1 \ d n+ \ 



Vn+l(«)> 



and the lowering-type formula: 



L (s,n)(p n (s) =-f n j^--j^(p n - 1 (s), 



(28) 



(29) 



where the raising-type operator L + (s,n) and the lowering-type operator L (s,n) are 
given by 



L + (s,n) = 



A„ r n (s) a(s) 



[n] q T' n Vx(s) 



1 



/+>/e( ._ 1W . ) iT, 



(30) 



and 



L (s, n) = 



[n], t£ [2n], Ax(s) 

+Ve(-M a + i)^^, 



(31) 



respectively. In the above formulas E~ f(s) = f(s — 1), E + f(s) = f(s + 1) and I is the 
identity operator. 

Notice that the last two formulas have a remarkable property of giving all the 
solutions (p n (s). In fact, from (31) setting n = and taking into account that (f-i(s) = 
we can obtain ipo(s). Then, substituting the obtained function in (30), we can find 
all the functions tpi(s),. . . , <p n (s),. . . . 

Proposition 3.1 The raising and lowering operators (30) and (31) are mutually ad- 
joint. 

Proof: The proof is straightforward. In fact using the boundary condition and after 
some calculations we obtain, in the case of discrete orthogonality, the expression 



6-1 



^2 ^n+lisi) 
Si=a 
6-1 



[2n] q 



A 



2n 



L + (si,n)cp n (si 



[2n + 2] ? 

L (si,n + l)ip n+1 {si) 



2n+2 



Ax( Si - ±) 

(p n (si)Ax(si - \) = a r 



dn+l 



The other cases can be done in an analogous way. ■ 

Proposition 3.2 The operator corresponding to the eigenvalue X n in (26) is self ad- 
joint. 

Proof: Again we will prove the result in the case of discrete orthogonality. Using the 
orthogonality conditions a(a)p(a) = a(b)p(b) = (which is a consequence of (13)), we 
can write 

b- 1 1 

<Pn(Si)y/e(Si ~ l)0-(s i ) Vx ^ ip l (s i - l)Ax( Si - ±) 

Si = d 

6-2 j 



s'.=o-l 
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^2 Vn(Si + 1) ^®(Si)(j(Si TTj — 1 Lpi(Si)Ax(Si + 5) + 

yn(Q)\/Q( Q - ^^^ v^aj ^'^ ~ l)Ax(a - §)- 



where in the last two sums we first take the operations A and V, and then substitute 
the corresponding value: e.g. Ax(a) = x(a + 1) — x(a). 

Now, we use the fact that (p n (s) = \J p(s)/d^ l P n (x(s)) q , as well as the boundary 
conditions a(a)p(a) = a(b)p(b) = 0, so 

y/e(a-l)a(a) ^ n {a)^i{a - 1) = ^/Q{b - l)u{b) Vn {b) V i{b - 1) = 0. 

The other terms can be transformed in a similar way. All these yield the expression 

b- 1 f 1 

^2 { V / Q( s ») cr ( g i + !) A / ^ n(Si + l)Ax(Si + |)+ 

b- 1 ( 1 

= ^2 Vn{Si) iy\/Q{Si)a(Si + l ) ^ x ( s .^ l( S i + l)Ax(Si + |) + 

^Q(si - l)a(si) ^ <fi(si - l)Ax(si - ±) j , 
from where the proposition easily follows. ■ 

4 Factorization of difference equation of hyper geometric 
type on the nonuniform lattice 

We will define from (26) the following operator 

H(s,n) ee Ve( S -lM*)^y£- + y/e{8)<T{8 + l)-^E+- 



O(s) , a(s) 



-X n Ax(s-^ I. 



Ax(s) Vx(s) 
Clearly, the orthonormal functions satisfy 

H(s,n)ip n (s) = 0. 

Let us rewrite the raising and lowering operators in the following way 

1 



L+(8, n) = u(s, n)I + V / 6MM«) y^)^' 



1 



L-( s , n ) = j)( s , n )/ + Ve(sW« + i)^£ + , 
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where, as before, 0(s) = a(s) + t(s)Ax(s — 5), and 

A n T„(s) O-(s) 



u(s, n) 



Ng t« Vx(s) : 



Proposition 4.1 T/ie functions u(s,n) andv(s,n) satisfy u(s + l,n) = v(s,n + l), or, 
equivalently u(s + 1, n — 1) = u (s, n) . 

The proof of the above proposition is straightforward but cumbersome. We will include 
it in appendix A. If we now calculate 

L~(s,n + l)L + (s,n) = v(s, n + l)u(s, n) + @{s)a{s + 1) ( A j( s ) ) + 

+u(s + l,n) |ve( S -lM S )^yi?- + V@(^ + l)^y^ + } , 

and substitute the values for u(s,n), v(s,n) and H(s,n) we get 

L~(s,n + l)L+(s,n) = /i T (n)I + u(s + 1, n)H(s, n), 
where the function 

fc=F(n) = f ilL^±l) _ ^±*U f - X n Ax(s - I)) + 

V[n], r' n Vx(s + l)J \[n] q r' n 2 V 

A TO T n ( S + 1) 9(g) 

[n] q T' n Ax(s) ' 

is independent of s. In fact, applying the last equality to the orthonormal function 
(f n (s) and taking into account (28) and (29), 



Similarly, 
where 



,qz, \ A2n A2n+2 

h{n) = sra;^ 1. 



L + (s,n — 1)L (s,n) = h (n)I + u(s,n — l)H(s,n), 



h±in) = ("S^ + - " - A '» + A " Al(s - i: 



[2n]/ 1 ' ™' Vx(s) 
K r n (s) \ 2n . . a \ ( 6(s - 1) 



[n] q r' n [2n\ q K > n >) \Ax{s-l) 

is independent of s. Furthermore, applying the last expression to the functions (p n (s), 
and taking into account (28) and (29), we obtain 

,+/ \ A2n-2 A2n 

n In) = 7 7-7 — 7- a„_i7„. 

V ; [2n-2] q [2n] q )n 

Remark: Notice that h^{n + 1) = h T (n). 

All the above results lead us to our main theorem: 
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Theorem 4.1 The operator H(s,n), corresponding to the hyper geometric difference 
equation for orthonormal functions ip n (s), admits the following factorization -usually 
called the Inf eld- Hull-type factorization- 

u(s + l,n)H(s,n) = IT (s,n + 1)L + (s,n) - h T (n)I, (32) 

and 

u{s, n)H(s, n + 1) = L + (s, n)L~(s, n + 1) - h T (n)I, (33) 

respectively. 

Remark: Substituting in the above formulas the expression x(s) = s we obtain the 
corresponding results for the uniform lattice cases (Hahn, Kravchuk, Meixner and Char- 
lier), considered before by several authors, see e.g. [9, 22, 15] and by taking appropriate 
limits (see e.g. [21, 25]), we can recover the classical continuous case (Jacobi, Laguerre 
and Hermite). 



5 Applications to some (/-normalized ortogonal functions 

For the sake of completeness we will apply the above results to several families of 
orthogonal (/-polynomials and their corresponding orthonormal g-functions that are of 
interest and appear in several branches of mathematical physics. They are the Al- 
Salam & Carlitz polynomials I and II, the big g-Jacobi polynomials, the dual q-Hahn 
polynomials, the continuous g-Hermite and the celebrated g-Askey-Wilson polynomials. 

The main data for these polynomials are taken from the nice survey [21] except the 
case of dual g-Hahn polynomials [3]. Nevertheless, they can be obtained also from the 
general formulas given in Section 2. 

Finally, let us point out that similar factorization formulas were obtained by other 
authors, e.g. Miller in [19] considered the polynomials on the linear exponential lattice 
and Bangerezako studied the Askey-Wilson case. Our main aim in this section is to 
show how our general formulas lead, in a very easy way, to the needed factorization 
formulas of several families for normalized functions -not polynomials-. 



5.1 The Al-Salam & Carlitz functions I and II 

The Al-Salam & Carlitz polynomials I (and II) appear in certain models of g-harmonic 
oscillator , see e.g. [4, 7, 8, 24]. They are polynomials on the exponential lattice 
x(s) = q s = x, defined [21] by 



t4 a) (^) = (-«)V 2 



- n ,x- 1 



291 



qx 
a 





and constitute an orthogonal family with the orthogonality relation (15) 



f Ui a \x; q)U$(x; q)u(x)d q x = d 2 J nm , 

J a 



where 

u(x) = (qx^^qxiq)^, and d 2 n = {-a) n {\ - q){q;q) n {q,a,a' l q;q) 00 q' 
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As usual, (ai, ■ ■ ■ , a p ; g) n = (a 1 ;q) n ■ ■ ■ (a p ; q)„, and (a; g), 
They satisfy a difference equation of the form (1) where 



„i/2 1 + a 

a ( x ) = (x - l)(x - a), r(x) = t(s) = r'x + r(0), being r' = f — -, r(0) = q 1/2 - 



1-q' 

The eigenvalues \ n and the coefficients of the TTRR are given by 



q-1 



A„, = \n] 



l-n/2 



9-1 

respectively. In this case we have 



and a n = 1, /3 n = (1 + a)g n , 7n = ag n x (g n - 1), 



a" = l, 5'(0) = -^±l 5(0) = a, <= y _ q 
The corresponding normalized functions (25) are 



2-" hnflll 
--, t„(0) = g 2 — -. 



<Pn(x) 



, (gx,a 1 qx-q) 00 (-a) n q\ 2 ) 
\ (1 - q)(q;q) n (q,a,q/a;q) 00 2 M 



gx 



Defining now the Hamiltonian for these functions (p n ( x ) 



^/a(x- 



and using that u(x, n) 



x 1 , v(x, n) = u(qx, n — 1) 



1-q 



1-q 



x 1 , thus 



L+(x,n) =u{x,n)I + q^- X r ^ ^ E~ , where £T/(z) = /(g _1 x), 



and 

L~(x, n) = v(x, n)I + 
we have 



x(q — 1) 
\fa(qx - l)(qx 



L (x, n + l)L + (x, n) = 



x(q - 1) 

aq l-n^ q n+l _ y 



E + , where E + f(x) = f(qx), 



and 



L + (x,n- 1)L~ 



\x,n) 



(q-lf 
aq 2 ~ n {q n - 1) 



I + v(x, n + l)H{x, n), 



I + u(x, n — l)H(x, n), 



(q-1) 2 

which give the factorization formulas for the Al-Salam & Carlitz functions I. If we now 
taking into account that (see [21, p. 115]) 

V^ a \x;q) = Ut ) (x;q- 1 ), 

then, the factorization for the Al-Salam &i Carlitz functions II 



(!) a*+ w ag;g oog( 2 ) 
\ (q,aq;q) s {l - q){q;q) n 



q ,x 



q" 

q; — 

a 



follows from the factorization for the Al-Salam & Carlitz functions I simply by changing 
q to g -1 . 
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5.2 The big g-Jacobi functions 

Now we will consider the most general family of g-polynomials on the exponential 
lattice, the so-called big g-Jacobi polynomials, that appear in the representation theory 
of the quantum algebras [31]. They were introduced by Hahn in 1949 and are defined 
[21] by 



D , 7 \ (aq;q)n{cq;q)n , [ q n ,abq n+1 ,x 

P n ( X ;a,b,c;q)= {ahqn+l . q)n 3^ agcq 

They constitute an orthogonal family, i.e., 

rag 



q;q , x(s) = q s = x. 



i cq 

where 

uj(x) 



i-aq 

I uj(x)P n (x; a,b, c; q)P n (x; a,b, c; q)d q x = d 2 n 8 n 

J CQ 

(a^x; q) (c 1 x;q) oo 



(x;g)oo(&c 1 x;q) 

oo 



d 2 = Qg(l ~ q){q,c/a,aq/c,abq 2 ;q) 00 (1 - abq)(q,bq, abq/c; q) n {-ac) n q 

n (aq,bq,cq, abq/c; q)oo (abq,abq n+1 ,abq n+1 ) n 

They satisfy the difference equation (1) with 

a(x) = q~ l (x — aq)(x — cq) and t(x) = t(x) = t'x + r(0), 

where 

, l-abq 2 1/2 a(bq - 1) + c(aq - 1) 
T = T, TT7? and T (°) = q ' 1 . 



and 

K _ _,-«/ 2[ „ ]s l_i^. 

I — q 

They satisfy a TTRR, whose coefficients are 
where 



(1 - aq n+1 )(l - cq n+1 )(l - abq n+1 ) 
(1 - abq 2n+1 )(l - abq 2n + 2 ) ' 

(1 -q n ){l -bq n ){\ -abc^q 11 ) 



q (1 - abq 2n ){l - a6(7 2n+1 ) 

Also, we have 

1 + abq 2 ~/ /r .N abq + acq + a + c 



a = 



-, <7(0) = , CT(0) = ac^, 



. q- n -abq n+2 tn . i-n a(bq 1+n - 1) + c(aq 1+n - 1) 
, T"n(0) '/ ■ 



n q y\i- q ) ' nv ; * 1-g 

The normalized big g-Jacobi functions are defined by 



/ (x/a,x/c;q) 00 (aq,bq,abq/c;q) 00 (abq,aq,aq,cq,cq;q) n (-ac) n ^ 
V (x, bx/c, c/ a, aq/c, abq 2 ;q) 00 (l - q)aq(l - abq) (q, bq, abq/ c; q) n 



392 



-n r,hr,n+l x 



aq, cq 



q;q 
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The corresponding Hamiltonian is 



yja(x-q)(x-aq){x-cq){bx-cq) y/ ' a{x-l)(x-a){x-c){bx-c) , 

H{x,n) = — hj +q — hj + 

x(q — l) x(q — l) 

l + abq 2n+l q(a + ab + c + ac) acq(q+l) _,\ . 

x 1 x I. 



q n (l — q) 1 — q 1 — q 



Furthermore, 



abq n+L acq 2 _ x afe<f t+1 acg _ x 

u(x,n) = x + -L> n x ana v(x,n) = x + D n _\ x , 

1 — q q — 1 1 — 9 <? — 1 

where 

_ a6(afe + ac + a + c)g 2n+3 - a(6 + c + a& + bc)q n+2 
n ~ (1 - abq 2n + 2 ){l - q) ' 

thus 

■ . , . . A/a(x — a) (x — aq) (x — cq)(bx — cq) ^ , . . „ . . 

L+(x,n) = u(x,n)I+-^ — — — —E~, where E~ /(x) = /(<T x), 

x(g — 1J 

and 

L-(x,n) =v(x,n)I + q-*—± A _ — where E + f (x) = f (qx) , 

x[q — i) 

so 

L~(x,n+ l)L + (x,n) = 5 n+ ij n+ il + v(x, n + l)H(x,n), 
L + (x, n — l)L _ (x, n) = 5 n ~/ n I + u(x, n — l)H(x, n), 

where 

_ (1 -a6g 2n " 1 )(l -a6g 2n+1 ) 



g2n-l(g-l)2 

The above formulas are the factorization formulas for the family of the big g-Jacobi 
normalized functions. 

Since all discrete g-polynomials on the exponential lattice x(s) = c\q s + C3 — the so 
called, g-Hahn class — can be obtained from the big g-Jacobi polynomials by a certain 
limit process (see e.g. [2, 21], then from the above formulas we can obtain the factor- 
ization formulas for the all other cases in the g-Hahn tableau. Of special interest are 
the g-Hahn polynomials and the big g-Laguerre polynomials, which are particular cases 
of the big g-Jacobi polynomials when c = q~ N_1 , N = 1, 2, . . ., and c = 0, respectively. 

5.3 The g-dual-Hahn functions 

In this section we will deal with the g-dual-Hahn polynomials, introduced in [3, 26] 
and closely related with the Clebsh-Gordon coefficients of the g-algebras SU q (2) and 
SU q (l, 1) [3]. They are defined on the lattice x(s) = [s] q [s + l] q by 



wCf ( , u\ _ (-l) n (g a - 6+1 ;g)n(g a+c+1 ;g)n , ( q- n ,q a - s ,q a+s+1 
w n {x(s),a,o )q - qn/2(3a _ b+c+1+n)k ^ q _ q ^ 3<p2y qa _ b+1 ^ a+c+1 

and satisfy a discrete orthogonality (14) with respect to the weight function 

q \((b-l) 2 -(2s-l)(o+c)) fgs-a+l ^ gS-c+1 ^ q s+b+l ^ g 6-s. q \ 



Pis) 



(1 - g)2(a+c-6)+l ( g> g> q s+a+l ? q s+c+l . ^ 
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where — \ < a <b — 1, |c| < a + 1, and for this weight function the norm is 

r ,U-4ab-4bc+6a+6c-8b+6+4n(a+c-2b)-n 2 + 17n+2b 2 ) /J>-c-n „\ 



(1 - g)2(a+ C -6+l)+3n [^(^ ^a+c+n+1 . ^ ' 

These polynomials satisfy a TTRR (17) with 

«n = l, 

/3 n =^( 2 ^+ 1 )[6-a-n+l] g [a+c+n+l] 9 +g5( 2 ^ 2 ^-^ 1 )[n] (? [6- C -n] (; + [a] g [a+l] 9 , 
7n = g ,2n+c+a_6 [a + c + n] 9 [6 - a - n] g [6 - c - n] g [n] g , 

1 n 

and the second order difference equation (1), whose eigenvalues are X n = [n] q q2~2 and 

a(s) = q^ s+c+a - b+ ^[s-a] q [s + b] q [s-c] q and t(x) = t(x) = r'x + t(0), 
with r' = -1 andr(O) = ^^ a - 6+c+1 )[a + l] 9 [6-c-l] g + g^ c - 6+1 )[6] 9 [c] g . 
Also we will need the values 
a" = k q , a'(0) = 2^(2[2], - qh~ b - - _ 

5 : (0) = 43(V^ + ^ 1 +^+2^ f> +2^ 1 - k ^-(l + ^)(^- 6 +g a + g c +g 1 - k ^^)), 
2re ? 

< = -<T n , r n (0) = ^( c - b -" +1 )[c+§] ? [6-§],+^(^ 1 -f)[a+§ + l] (? [6-c-n-l] ? , 

In this case, the Hamiltonian, associated with the g-dual Hahn normalized functions 

y / p(s)/dlW^(x(s);a,b) q , is 

1{ + ^ ♦ /([* + l] 2 " [a]2)([6]2 - [s + 1]2)([ S + 1]2 _ [c] 2) 

#(s,n) = ^+^+2) V ? 1 J! « « 

[2s + 2] g 



^+^+ 2 ) V £?- - g*"* [n],[2s + 1],/+ 

|( c+ a-6+2) ( [s ~ a}g[s + b} q + [s - c] q [s + 1 - a] g [s + 1 + + 1 - c] q \ 
Q \ [2a], [2 S + 2] q ; ' 

where = /(s + 1) and £T/(s) = /(s - 1). Then, using that 

u(s,n) = q*-$x(8 + n/2) - q\+^ ( q ¥ c - b ~ n+1 )[ c + § ],[& - 

^(a+c-b+l-f + n + 1]?[6 _ c _ n _ l]q) _ q ± {s+ c + a- b+ 2) [s - a] q [s + b] q [s - c] q ^ 

[2sjg 

and taking into account that v(s, n) = u(s + 1, n — 1), we find 



, V '(M2 - W|)([6]i - Wi)(Wl - [c\D 

[2sJ<j 



and 

1{ + ^ */([* + I] 2 - M 2 )([6] 2 - [a + 1] 2 )([, + l] 2 - [c] 2 ) 
L-( S , n) = v(s, n)I + ^(^+ 2 ) V ___ /, . 
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Thus 

L~(s, n + l)L + (s, n) = q~ 2n j n+1 I + v(s, n + l)H(s, n), 

and 

L + (s, n - l)L~{s, n) = q- 2n+2 j n I + u(s, n - l)H(s, n), 
are the factorization formulas for the g-dual Hahn normalized functions. 

5.4 The Askey Wilson functions 

Finally we will consider the family of Askey- Wilson polynomials. They are polynomials 
hi 



on the lattice x(s) = ^(q s + q s ) = x, defined by [21] 



(ab;q) n (ac;q) n (ad;q) n ( q n ,q n 1 abcd,ae l9 ,ae l 

p n {x(s);a,b,c,d) = 4 3 

a™ \ ab, ac, ad 



q;q 



i.e., they correspond to the general case (23) when q Sl = a, q S2 = b, q ss = c, q S4 = d. 
Their orthogonality relation is of the form 

/i 
oj(x)p n (x;a,b,c,d)p m (x;a,b,c,d)y/l — x 2 K q dx = Snmd^, q s = e ld , x = cos9, 



where 



2nn q {l — x )n{x,a)n{x,b)n(x,c)n{x,a) 
and the norm is given by 

(abcdq n ~ 1 ;q) n (abcdq 2n ; q)^ 



dt = 



(q n+1 ,abq n , acq n , adq n , bcq n , bdq n , cdq n ; q)^ 



The Askey- Wilson polynomials satisfy the difference equation (1) with 

<*{s) = -q- 2s+1/2 K 2 q (q S ~ a)(q s - b)(q s - c)(q s - d), n q = (g* - q ~h) 

and t(x) = t(x) = r'x + r(0), where 

T ' = 4(q - 1)(1 - abed), r(0) = 2(1 - q)(a + 6 + c + d - abc - abd - acd - bed). 

Furthermore, they satisfy the TTRR (17) with coefficients 

R a + a~ l - {An + C n ) _ CnAn-l 

a n — A) Pn — 2 ' T n — ^ ' 

where A n , C n are defined by 

(1 - c% n )(l - acg n )(l - a^ n )(l - abedq 11 " 1 ) 
n ~ a(l - a6cdg 2n - 1 )(l - abcdq 2n ) ' 

_ a(l - g ra )(l - bcg"- 1 )(l - M^Xl - edg^ 1 ) 
n ~ (1 - abcdq 2n - 2 )(l - abedq 2 ™- 1 ) ' 
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and whose eigenvalues are X n = 4q n+1 (l — q n )(l — abcdq n x ). In addition, we have 
a" = -4(q - 1) 2 (1 + abcd)q~ 1 ' 2 , 

a'(0) = (q- l) 2 (a + b + c + d + abc + abd + acd + bcd)q~ 1/2 , 
5(0) = (q ~ 1) 2 (1 — ab — ac — ad — be — bd — cd + abcd)q~ 1 ^ 2 , 
T^ = 4q- n (q-l)(l-abcdq 2n ), 

t„(0) = 2{q-l){-a-b-c-d+{abc + abd+acd + bcd)q n )q- n/2 . 



Defining now the normalized functions (see (15)) y / co(x)/d 2 l p n (x; a, b, c, d), the corre- 
sponding Hamiltonian H(s, n) is 

2 3/2 2o 3 / 2 
H(s, n) = p g _y G{s, a, b, c, d)E~ + ^—^-G{s + 1, a, b, c, d)E + + 

2 1 - 2s+ i/ 2 nti(i-^ +s ) , - 2s+V 2 ny^o , 

^ [2a + 1], + 9 [2a -1], 

g-" +1 K 2 (l - g n )(l - abcdq n - l )[2s]^jl 

where 



G(a, a, b, c, d) 



H(l - 2q°*q-V*x{8 - 1/2) + q~ l q 2 ^), 



We now define 



where 

D n = -Aq' n l 2+1 ' 2 {q - 1)(1 - ated^ 1 ). 

_ (-a - b - c - d + (abc + abd + acd + bcd)q n )q n/2 
n ~ 2(1 - abcdq 2n ) ' 

Taking into account that v(s, n) = u(s + 1, n — 1), we find 

2<7 3 / 2 

L (s, n) = u(s, n)I + — — G(s,a,b,c,d)E , 

[2s — ljg 

2<? 3 / 2 

L (a,n) = v(s,n)I + j- G(a + 1, a, 6, c, gQ£ + , 

[za + l\q 

where £ _ /(a) = /(a - 1) and E+f(s) = /(a + 1). Thus, 

L~ (a, n + l)L + (a, n) = D 2n D 2n+ 2ln+iI + v(s, n + l)H(s, n), 

and 

L + (s,n - l)L~(s,n) = D 2n - 2 D 2n ^ n I + u(s, n - l)H(s,n), 
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which is the factorization formula for the Askey-Wilson functions. 

To conclude this paper let us consider the special case of Askey-Wilson polynomials 

when a = b = c = d = 0, i.e., the continuous g-Hermite polynomials 



H n (x\q) = 



\0 



q; q n e~ 2ie 



X = cos ( 



These polynomials are closely related with the g-harmonic oscilator model introduced 
by Biedenharn [13] and Macfarlane [23], as it was pointed out in [8], where the factor- 
ization for the continuous g-Hermite polynomials were considered first. If we substitute 
a = b = c = d = 0in the above formulas, we obtain the factorization for the g-Hermite 
functions 



! h(x, l)h(x, -l)h(x, q^hix, -q 1 / 2 )(q n+1 ;q) c 



I 2lTK q (l-X 2 ) 

In fact, since for continuous g-Hermite polynomials 



■H n (x\q). 



a(s) 



2 2s+l/2 



-K,q 



t(s) = % - l)x(s), \ n =4q~ 



-n+l. 



1 



and the coefficients for the three-term recurrence relation are a n = 1, fi n = 0, 7„ 
(1 — q n )/4:, then we obtain 



2o 3 / 2 



2o 3/2 f -28+1/2 2s+l/2 
-E + + 2\^ — — + 



[2s + l\ q 



[2s + l] q [2s - l] q 



-q- n+1 K 2 q (l-q n )[2s} q )l, 



L+(s,n) = \ -4q~ n / 2+1 / 2 (q - l)x(s + n/2) + t^'l 1 /+ \ E 



[2s - 1], 



2s+5/2 



L-{ 8 , n) = ^-4q- n / 2+1 (q - l)x(s + n/2 + 1/2) + J-^ 
and h ± (n) = 4K 2 q- 2n+1 (l - q n ). 



[2s - 1] 5 

[2S + l]g 
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Appendix A 

Here, for the sake of completeness, we will prove Proposition 4.1, by showing that u(s + 1, n) ■ 
v(s, n + 1) = 0. To do that, we start with computing the difference 



u(s + 1, n) - v{s, n + l) — — ; : — —+ 



n, 



X n +1 T n+ l(s) 

[n + l] 



A„ + iAa;(s - \ 



r' n Ax(s 

^2ra+2 



<? 'n+l 



[2n + 2] ( 



(x(s) - fin+l) + 



Ax(s) 
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Now we use the expansion r n (s + 1) = r' l x n (s + 1) + r„(0). Since 

A^) = a; 2( s + i)_ a; 2 ( ) = ^ + = C^?" 1 + 1) + 2C 3 = 

Aa;(s) a;(s + 1) - a;(s) 
(Cig s +* + C 2 g- s -*)[2], + 2C 3 = [2} gXl (s) + (2 - [2],)C 3 , 

x( S )Ax( S - 1) = a;(s)(C l9 -*(g - 1) + Cag-^^r?- 1 - 1)) = x( S )(C 1<Z s - C 2 <T s )fc 9 - 

(CjV 8 - C| <Z - 2s )fc 9 + C 3 (C iq s - C 2 q- S )k q , 
where k„ = q? — 



Ax(s 



(x(s)Ax(s-±)) 



'(C 2 q 2s+1 + Ctq- 2s - 1 )[2\ q + Cz(C 1 q s+1 2 +C 2 q- s - 1 i) 
Ciq s+1 i -C 2 q- s -i 



and 



Then 



Ax(s) 



(A^-^))~(W-C 2 ^) fcg ) 



C iq s+ i - C 2 q- s -i 



-kg. 



Aa(s) A 
Ax(s) ~ Ax(s) 



(a(s) \r(s)Ax( S ±)) = 



Ax(s) 



^x 2 (s) + a'(0)x(s) + 5(0) - \ (t'x(s) + r(0)) Ax(s - ±) ) = 



Y « 2 W*) + ( 2 - [2],)C 3 )+?'(0) - |r(0) ' ( V ^ + ' 



1 , / [2] q (C*q 2s+1 + C'h- 2 *- 1 ) + C 3 (C iq s +? + C 2 q- S -^) 



2 T 



C iq s+ 2 - C 2 q- s -2 
This yields for u(s + 1, n) — v(s, n + 1) the expression 



q 



^[2W S ) + ^(2-[2] g )5" + a'(0y 



r' ([2]q{Clq 



2 „2s+l 



■CZq 



2„-2s-l 



2 

r(0) / x 1 (s)-C 3 
2 VCi<z s+ ^ -C 2 g- S -^ 

A 2 n+2 



+ 



C sXl {s)-Cl 



C iq s +i - C 2 q~ s -2- C iq s +? - C 2 q- s -i) q 



X n +1 T n+1 (s) 

[n+1] 



A„+iAa;(s - |] 



q 'n+l 



[2n + 2] 9 
[n + 2] g r„+i(0) 



c lg - + c 2g - + c 3 - [n + 1] f" (0) + 



'n+l 



C 3 (l + [n + l], 



■2],) 



t( S )Ax( S ~1) 
Az(s) 



Next, we expand Ai„(s) and ^-[2]xi(s), make some straightforward calculations and use the 
identities: 



A n t„(0) A 2n+2 r„(0) 

M, < +[n + lk [2n + 2] g r' n 



An . r -.1 A 2 „ +2 

j— j- + [n + l] q -, 



[2n + 2] t _ 



T„(0) 



[n + 2],T„(0), 



A„+i r„+i(s) A 2n+2 r n+ i(0) A„+i 
— rTi 1 t n + 2 ]<?-^ — ; = [n + l] 9 T n+ i(0) + — — r x n+1 (s), 



iq 'n+l 



[2n + 2] 



9 'n+l 



1],' 
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as well as 



AL(C l9 s + 1+ t + C^— 1 "*) - ^-[2} q (C iq s+ i + C 2 q- S -?) - X2n + 2 (C iq s + C 2 q- S )- 

[Tt\q Z [Zfl "T" ^Jq 



Xn+l(C iq s - C 2 q- S )k q + ±r'(C l9 s+ * + C 2 q- s -^)(q + q' 1 ) = + q- 



C iq s a" . n+ i i C2<rV, _ n _i _± C 2 <r*?" _„_i _i 

( g «+ 2 + (q 2 +q 2 )+ -rr(-g 2 + g 2 ) = 



2(<p —52) 2 2(<p — q 2) 

[n + l]g 

we find 

= - r^rTT^i^ 2 * 1 + C 2 q- S -^) + C 3 -^ - [n+ 2] 9 r„(0) - C 3 a" - ?'(0) + \r'C z k q + 
[n + i\ q [n\ q z 

K(0)k q + {n+l} q T n+1 (0) + j^^^ 

Finally, we substitute the expression for r„(0) and use the identities 
-[n + 2} q [n] q - l + [n+l] q [n+ l] q = 0, 
-[n + 2} q (q n / 2 + q- n / 2 ) + k q + [n + l] q (q^ +1 ) ' 2 + q^+D/ 2 ) = 0, 
and the result follows. 
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